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Abstract 

A non-dissipative model for vortex motion in thin snpercondnctors is considered. The La- 
grangian is a Galilean invariant version of the Ginzburg-Landan model for time-dependent 
helds, with kinetic terms linear in the hrst time derivatives of the helds. It is shown how, 
for certain valnes of the coupling constants, the held dynamics can be reduced to hrst order 
diherential equations for the vortex positions. Two vortices circle around one another at 
constant speed and separation in this model. 
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1. Introduction 


Magnetic flux penetrates a Type II superconductor in the form of vortices [1], and recently 
it has become possible to produce images of vortices sufficiently rapidly that their mo¬ 
tion can be observed directly [2], In the Ginzburg-Landau theory of superconductivity, a 
charged scalar held representing the electron-pair condensate is coupled to the electromag¬ 
netic held. The basic vortex solution, discovered by Abrikosov [3], is a localised magnetic 
hux tube surrounded by a circulating supercurrent. 

The Ginzburg-Landau potential energy functional contains only one dimensionless cou¬ 
pling constant A. The value A = 1 (in our units) is mathematically particularly interesting, 
because in this case there are no forces between static vortices, and there is a continuous 
family of static multivortex solutions. A Type II superconductor is modelled by A > 1. 
In this case, the potential energy of a two-vortex conhguration decreases as the separation 
increases, in other words, vortices repel [4]. However, there are several possibilities for how 
the vortices might move, depending on the nature of the dynamical equations for the helds. 
Let us ignore pinning, which tends to prevent vortex motion at all. The hrst possibility 
is that the vortex acceleration is proportional to the force acting. This is what occurs 
in the relativistic generalisation of the Ginzburg-Landau model, known as the Abelian 
Higgs model. Relativistic vortices may be interpreted as a solitonic version of fundamental 
strings [5], or as strings joining conhned quarks, or as cosmic strings produced at a phase 
transition early in the universe’s history [6]. The second possibility is that the vortex 
velocity is proportional to the force. This is modelled by dissipative equations involving 
the hrst time derivatives of the helds [7]. Recently, one version of such equations, the 
Ginzburg-Landau gradient how equations, have been analysed in detail [8]. The third 
possibility is that the vortex motion is at right angles to the force, so two vortices circulate 
around each other, as in a huid, and there is no dissipation [9]. In real superconductors, 
vortex motion is usually dissipative, but at very low temperatures, it has been argued that 
the third type of motion would occur [10]. 

The response of vortices to an applied electric “transport” current, perpendicular to the 
vortex cores, can distinguish the second and third types of motion. If the vortices move at 
right angles to the current, then the dynamics is dissipative, but if they are carried along 
by the current (again, as in a fluid) then the dynamics is non-dissipative. To see this, note 
that a moving vortex has an electric held in its core, perpendicular to both the velocity 
and the direction of the magnetic hux. Also, part of the applied current penetrates the 
vortex core, where it becomes a normal electric current. There is dissipation when the 
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current and electric field in the core are parallel, but not when they are perpendicular. 

The purpose of this paper is to analyse a model for the held dynamics in a thin, essentially 
two-dimensional superconductor, and to show that it leads to vortex motion of the third 
type. The held equations are obtained from a Lagrangian whose kinetic terms are linear 
in the hrst time derivatives of the helds and whose potential part is the usual Ginzburg- 
Landau energy. This Lagrangian and its associated held equations are given in section 2. 
The Lagrangian is Galilean invariant, so we can see precisely how vortices respond to a 
transport current. 

Section 3 is a review of the static vortex solutions of the Ginzburg-Landau theory, focussing 
especially on the manifold the 2n-dimensional parameter space of static n-vortex 
solutions which exist at the critical coupling A = 1. These solutions are also present in the 
model considered here, at special values of the couplings. 

Section 4 treats the case where A is close to 1, and where n-vortex motion can be well 
approximated by a motion on M'^. This adiabatic approximation assumes that at each 
instant the held is a static solution, but that the parameters of the static solution (i.e. the 
vortex positions) slowly vary with time. It is shown how the kinetic energy and potential 
energy of the held Lagrangian can be reduced to give a hnite-dimensional dynamical system 
on M^. The kinetic energy can be expressed in terms of local data associated with each 
of the vortices, and although an explicit form for this is not obtained, some conclusions 
can be drawn. In section 5 the equations of vortex motion are derived from the reduced 
Lagrangian, and it is shown that two vortices circle around one other. 

The method used here for reducing the held dynamics to a particle dynamics for vortices 
closely follows the analysis of slow vortex motion in the Abelian Higgs model with A = 1, as 
carried out by Samols [11]. There the reduced system’s kinetic energy expression involves a 
Riemannian metric on the potential energy is a constant, and the vortex dynamics is 
modelled by geodesic motion on [12]. The geodesic motion is modihed by the ehect of 
a potential energy varying over M” if A is not exactly unity. This was studied by Shah [13]. 
Our reduced system has the same potential energy, and a kinetic energy which involves a 
connection or gauge potential on that depends on the same data as Samols’ metric. 

2. The Schrodinger—Chern—Simons Lagrangian 

Let us consider a two-dimensional superconductor, with translational symmetry. A thin 
him with no defects might be close to this idealization. There is a complex scalar held 
(/), representing the electron-pair condensate, and a gauge potential with time component 


3 



at and spatial components at : i = 1,2. The magnetic field is B = dia 2 — d 2 ai and 
the electric field Ei = dtat — di- (An overdot denotes ^.) Units are fixed so that (f) has 
magnitnde 1 in the condensed phase (vacnnm), and its covariant time and space derivatives 
are Dt(t> — dt4> — and Di(f) — di(f) — iatcf). 

The Lagrangian that we shall consider, L, is an expression involving no higher than the 
first power of time derivatives. Explicitly, 


L 



+ n {Bat + E 2 ai - Eia 2 ) - jat 



- j(l - li-l")" - aiJj)d^x , 


( 2 . 1 ) 


where \B^ + ^{Dicj))*Dicj) + -1(1 — is the standard Ginzbnrg-Landan field energy 

density for static fields. Since there is no relativistic invariance, the snmmation convention 
is nsed only in the two space dimensions. The term with coefficient // is the Chern- 
Simons density for the gange field. 7 , fi and A are real constants, with A positive. Jj 
is the transport cnrrent, that is, a constant vector in the plane of the snpercondnctor. 
This Lagrangian is hardly original. The scalar field part of the Lagrangian conpled to the 
Chern-Simons term has appeared in the theory of Chern-Simons vortices, developed by 
Jackiw, Pi and others [14]. The inclnsion of the yat term was advocated by Barashenkov 
and Harin to allow the possibility of a condensate (|(/)| = 1) at infinity [15]. Somewhat 
original is that the term is inclnded, bnt not the ^EtEt term nsnally present in the 
Maxwell Lagrangian. This is jnstified becanse there is no relativistic invariance in the rest 
of the Lagrangian, and it is desired to have only first time derivatives in the field eqnations. 

The field eqnations are 

i^Dtcj) = ~DiDi(j) - ^{1 - \(j)\y(j) ( 2 . 2 ) 

eijdjB = jf - jy + 2neijEj (2.3) 

2 /iS = 7(1 - |(/>|2) , (2.4) 


obtained by varying with respect to (j)*, at and at respectively. Jf is the snpercnrrent, 
defined by 


jf = -'-{rDA-'i>(D,4,r). 


(2.5) 


Eq.(2.2) is the ganged non-linear Schrodinger eqnation, eq.(2.3) is Ampere’s law, a two- 
dimensional version of V x B = J, and eq.(2.4) is a constraint. There are really two 
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contributions to the total current, namely, the super current Jf and the normal current 
j/^ which is present only in the cores of the vortices. This being a non-dissipative model, 
there is no Ohmic conductivity. However, there is a Hall conductivity 2//, and hence = 
2iieijEj. The effect of a Hall conductivity on vortex motion was previously considered by 
Dorsey [16], and equations (2.2)-(2.4) are just a special case of those analysed by Dorsey. 

It is possible to have a constant, asymptotic value for the supercurrent jf. This asymptotic 
current is called the transport current Jj. An everywhere constant current arises if = 
exp i(k ■ X — lkpt/ 27 ), for example, and then = k. Such a constant current in a thin 
him produces a magnetic held above and below the him (^k x z above and — ^k x z below) 
but no held in the him, and the three-dimensional Ampere’s law is thereby satished. The 
version of Ampere’s law (2.3) leaves out d jdxj, terms. To correct for this and avoid a linear 
growth of B in the plane, due to the constant current, it is necessary to subtract oh 
from J*^, as in (2.3). This also explains the need for the —aiJf term in the Lagrangian 
density. 

The non-linear Schrodinger equation and its complex conjugate imply the supercurrent 
conservation law 

■ ( 2 - 6 ) 

In addition, there is Faraday’s law 


diE2 — d2Ei + — 0 (2.7) 

which is an immediate consequence of the dehnitions of Ei and B in terms of the gauge 
potential. (The equation V ■ B = 0 is trivially satished in the two-dimensional geometry.) 
Taking the divergence of (2.3), and combining (2.6) and (2.7), one hnds 

k(2MB + 7H")=0. (2.8) 

To avoid explicitly breaking translational invariance, 2iiB -\- 7 |(/>P must be a constant, 
independent of position; to admit the condensate |(/)| = I with no magnetic held, this 
constant must be 7 . Hence 2ij,B = 7(1 — |(/>P), which agrees with (2.4). 

An interpretation of eq.(2.4) is that the total electric charge density is zero. The charge 
density due to the condensate is — |(/>P, and the charge density of the background positive 
ions is I. Where the condensate is absent, the background positive charge is neutralised by 
decoherent, normal electrons; the normal charge density is —^B. However, this is only 
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an approximation, and the electric charge density is not exactly zero. Note that there is 
no eqnation for diEi among onr held eqnations; this is becanse of the absence of ^EiEi in 
the Lagrangian. One may evalnate the charge density p nsing diEi = p; it is expected to 
be very small. Since Ampere’s law reqnires that = 0 , one cannot nse the cnrrent 

conservation eqnation of electrodynamics to dednce anything abont p. Physically, there 
will be small electric charges generated by moving vortices, and diJl°^^^ will not be exactly 
zero, bnt this can be ignored in the non-relativistic limit. 

A remarkable property of the system of eqnations (2.2)-(2.4), together with (2.6) and 
(2.7), is that they are exactly Galilean invariant. It was stressed by Aitchison et al. 
that the eqnation for the scalar held shonld be Galilean invariant [10]. Onr additional 
eqnations dehne a non-relativistic limit of Maxwell’s eqnations with the same invariance. 
In fact, we have a “magnetic version” of Galilean invariance, with snbstantial cnrrents 
and negligible electric charge density. For an illnminating discnssion of Galilean invariant 
limits of electromagnetism, see ref. [17]. 

The basic Galilean transformation is a boost by a velocity v. Gange invariant scalar 

qnantities transform as / —;► /' where /'(x, t) = /(x — vf, t). Let ns denote x' = x — vt. 

The held transformations are 

(/)'(x, t) = (/)(x', ( 2 . 9 ) 

a((x, t) = ai{x.\t) ( 2 - 10 ) 

a((x, t) = at(x', t) - niai(x', t) . ( 2 . 11 ) 

It is straightforward to verify that eq.(2.2) is invariant nnder these transformations, in the 
same way that the Schrodinger equation in quantum mechanics is Galilean invariant when 
there is no potential. The magnetic and electric helds transform to 


S'(x,t) = 5(x',t) (2.12) 

F;'(x, t) = Ei{-x' , t) - CijVjBi^x', t) (2.13) , 

and the supercurrent becomes 

J'^(x, t) = J^(x', t) + 7 v|(/>(x , t)\^ . (2.14) 

The transport current has therefore also transformed to 

J''^ = J^ + 7v. (2.15) 
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Combining eqs.(2.12)-(2.15), we find that eq.(2.3) is Galilean invariant provided the con¬ 
straint (2.4) is satisfied. Finally, it is easily checked that the constraint (2.4) itself remains 
unchanged by a Galilean transformation. 

The physical interpretation of the Galilean invariance is as follows. Given any solution of 
eqs.(2.2)-(2.4) in the absence of a transport current, the effect of a transport current is 
simply to boost the solution so that it drifts along with the current at a velocity v = ^J^- 

Since we now understand the effect of a transport current, let us from now on assume there 
isn’t one, and that the fields have finite energy, with all currents and electromagnetic fields 
localized in a finite region of space. 

3. Vortices 

In the Ginzburg-Landau theory, fields which are smooth and of finite energy have the 
asymptotic behaviour \(^\ 1 and Dicf) —0 at spatial infinity. Such fields are classified 

by their integer winding number. When the winding number is n, the phase of cf) increases 
by 27rn anticlockwise around the circle at infinity. The vanishing of Dicf) at infinity implies 
that the gauge potential also carries information about the winding number, from which 
follows the magnetic flux quantization 

f B (fx — 27rn . (3-1) 


If the zeros of (f) are isolated, then the net number of zeros, counted with their multiplicity, 
is also n. A zero of multiplicity 1 may be identified as a magnetic flux vortex, and a zero 
of multiplicity —1 as an antivortex, so the winding number is also the net vortex number. 

The static Ginzburg-Landau equations 

AA</>+^(l-!</>!")</> = 0 (3.2) 

eijdjB = jf (3.3) 


have a vortex solution with unit winding number of the form (in polar coordinates) 


Cli 


(f) — k{r)e^^ 


r 


sin 6* , 


0-2 


r 


cos 9 , 


(3.4) 


where k ^ 1 and <7 —^ 1 exponentially fast as r ^ oo, and where k and g both vanish at 
r = 0. The precise form of k and g must be determined numerically. More generally, for 
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winding nnmber n > 1 , there are mnltivortex static solntions of a similar form, bnt with 
(f) = /c(r)e*"^^, and where /c —1 and ( 7 —s>nasr—>cx). IfA>l, vortices repel, so the 
circnlarly symmetric mnltivortex solntion is nnstable to break-np into individnal vortices. 

In the special case A = 1, there are no static forces between vortices. As Bogomolny showed 
[18], in this case eqs.(3.2) and (3.3) are satished, and the energy minimised, provided the 
hrst order Bogomolny eqnations 

Di(f) + iD2(j) = 0 (3.5) 

(3.6) 

are satished. Eq.(3.5) only permits cf) to have zeros of positive mnltiplicity, so solntions 
consist only of vortices. (To get antivortices, the sign in (3.5) and the sign of B shonld be 
reversed.) The Bogomolny eqnations have not only circnlarly symmetric solntions, bnt, for 
any n > 1, solntions with vortices located at n arbitrary points in the plane. More precisely, 
Tanbes proved that, modnlo gange transformations, there is a nniqne hnite energy solntion 
of eqs.(3.5) and (3.6) with (f) having zeros at any n prescribed points {x’' : 1 < r < n}, 
some of which may coincide [19]. Each snch solntion has winding n, total hnx 27rn, and 
energy im. 

The parameter space of snch solntions, which is called the n-vortex modnli space is 
a smooth manifold of dimension 2n. To see this, identify and C, and regard the zeros 
of (f) as the complex nnmbers {zr : 1 < r < n}. Geometrically, M'^ is the manifold C”/S” 
where is the permntation gronp acting on the n zeros of (f). /TB is actnally smooth 
despite the apparent conical singnlarities where two or more zeros of (f) coincide. The zeros 
nniqnely dehne a polynomial P{z) = z'^ + aiz'^~^ + ... + with precisely these zeros. 
The coefficients {ai,..., a^}, which are symmetric polynomials in the zeros, are arbitrary 
complex nnmbers, so, as a manifold, is also the space of monic (leading coefficient 
= 1) polynomials of degree n, and this is simply C"^. The coefficients {ai,..., a^} rather 
than the nnordered zeros {zi ,..., Zn} are the “good” coordinates on 

The eqnations (2.2)-(2.4) that we are interested in are not simply the static Ginzbnrg- 
Landan eqnations. However, if A = 1 and /U = 7, any solntion of the Bogomolny eqnations 
(3.5)-(3.6) is also a static solntion of eqs.(2.2)-(2.4). The constraint (2.4) is one of the 
Bogomolny eqnations, so it is satished. Dtcf) and Ei can be consistently set to zero, with 
at = 0 . 

One might seek static solntions of eqs.(2.2)-(2.4) for A 7 ^ 1. Presnmably snch solntions 
exist, if at is allowed to be non-zero. They represent stationary points of the Ginzbnrg- 



Landau energy, subject to the constraint (2.4), with at a Lagrange multiplier held. Solu¬ 
tions of the Bogomolny equations will come close to being solutions if A ~ 1 and ~ 7 . 
As in the earlier case, with A 7 ^ 1 , only circularly symmetric static solutions are expected, 
and they may again be unstable. 

The constraint (2.4) is very far from being satished by the static solutions of the Ginzburg- 
Landau equations in the extreme Type II regime (A 3> 1). It is therefore unclear whether 
our model is of any relevance to vortex dynamics in this regime. Instead, we shall consider 
the case of A ~ 1, which is realised by niobium and vanadium in certain temperature 
ranges. 

4. A Reduced Lagrangian 

The rest of this paper is devoted to constructing a reduced Lagrangian for n-vortex dy¬ 
namics, assuming A ~ 1. We need to assume that /U = 7 ; this is essential for simplifying 
the kinetic energy. The method is similar to that used to study n-vortex dynamics in the 
Abelian Higgs model, both at A = 1 and when A ~ 1. Let us consider helds which at each 
instant are static solutions of the Bogomolny equations, but where the moduli, that is, the 
vortex positions x’’ (or better, the coefficients of the polynomial P{z)) are time-dependent. 
These moduli will vary slowly if A is close to 1. The helds are inserted in the Lagrangian 
(2.1), and the integrals carried out, where possible. The result is a reduced Lagrangian for 
motion on the moduli space The reduced system is an approximation to the true held 
dynamics, but we shall not try to estimate the errors involved. The use of solutions of the 
Bogomolny equations is possibly better justihed here than in the context of the Abelian 
Higgs model, because the constraint (2.4), which is one of the Bogomolny equations, must 
be satished. 

Let us denote by {A" : 1 < a < 2n} some general coordinates on the moduli space 
for example, the components of the vortex positions. The held Lagrangian has a kinetic 
term which is hrst order in time derivatives, and a potential term. The reduced system is 
therefore expected to have a Lagrangian of the form 

L = A«(X)A"-I/(X) . (4.1) 

L may have additional terms which are total time derivatives, but these do not ahect the 
dynamics. has the interpretation of a gauge potential or connection on and it is 
somewhat arbitrary since replacing it by = Aot+daA (a gauge transformation) changes 

L by the total time derivative From L we hud the equations of motion 

+ dpV = 0 (4.2) 
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where Bct /3 = d^Ap — df^Aa is the curvature of the connection A. The motion according 
to (4.2) is non-dissipative, with V constant along any solution path 

The potential energy in the reduced model is the Ginzburg-Landau potential energy, but 
because the helds satisfy the Bogomolny equations, this simplihes to [18] 

V = mi + J cPx . (4.3) 

The constant mi does not affect the dynamics. The integral without the factor i(A — 1) 
is positive, and is some function of the relative positions of the vortices, invariant under 
a rigid rotation. The detailed form of the integral is not known, for general n, but it 
is expected to be minimal when the vortices are well separated, and maximal when the 
vortices are coincident. The integral has been computed in the case of two vortices, by 
Shah [13]. It increases monotonically as the vortex separation decreases. 

We can proceed much further with the calculation of the kinetic energy in the reduced 
model. Let us start with the kinetic terms of the held Lagrangian, with 7 = // 

T = 7 J ~ - at + Bat + {d 2 at - ^ 2)01 - {diat - di)a 2 ^ (fx . (4.4) 

It helps to make some assumptions about the asymptotic gauge. Recall that vortices are 
exponentially localized. Gauge invariant quantities B, {Dtipl and \4>\ approach their asymp¬ 
totic values exponentially fast. Suppose that the solutions of the Bogomolny equations are 
in the gauge where, for large r, they are of the form cf) — ai — — ^ sin^, a 2 = y cos 0 , 
with at most exponentially small corrections. Since Ei and Dtcf) are exponentially small 
asymptotically, in this gauge at is exponentially small too. 

Using the appropriate version of Stokes’ theorem, the terms involving derivatives of at can 
now be removed, and T expressed as 

T = -f j -|- a 2 di — aid 2 ^ d^x + 7 ^— 1 + 2B)at d^x . (4.5) 

The boundary terms vanish as a result of the gauge hxing, and the Bogomolny equation 
(3.6) (or the constraint (2.4)) implies that the second integral vanishes. Now dehne a = 
ai + ia 2 - Then T reduces to the compact expression 

T = -7 / + a*a) d^x . (4.6) 
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Let us next express the scalar field as 


^ ^ Q^h+ix (4 7) 

where h and x are real. (/> is a smooth function, but the gauge is not fixed except asymp¬ 
totically, so there is considerable freedom in the choice of x- Because of the presence of 
vortices, x is multivalued, increasing by an integer multiple of 2 ti along an anticlockwise 
loop around any zero of (f). However, in a given gauge, the gradient and time derivative of 
X are well-defined. 

Away from the zeros of (/>, the Bogomolny equation (3.5) implies that 

a = ^d2h + dix + ii-^dih + d2x) (4.8) 

and hence 

a = ^d2h + dix + + ^sx) 

Also 

Therefore 

T = -xj (e'^X + i^d 2 h + dix){-^dih + d 2 x)-i-^dih + d 2 x)i^d 2 h + dix)^ (fx. (4.11) 
It is easy to verify that (4.11) can be reexpressed as 

r = -7 /((e* - • Vx) + az/i - ai/ 2 ) d^x (4.12) 

where 

fi = {d2h + dix)x + ^dihh, f2 = {-dih + d2 x)x+^d2hh. (4.13) 
Since, from (4.8), B = —the Bogomolny equation (3.6) reduces to 

V^h-e’^ + 1 = 0 , (4.14) 


(4.9) 

(4.10) 


so T simplifies to 


T = -x 


1 O 

(x- 2 gj(VA-Vx) + 52/i-5i/2) d'‘x. 


(4.15) 
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To progress, we need to specify carefully the region of integration. It will simplify our 
calculations to assume that at any time t, the vortex positions x^(t) are distinct, so (f) has 
a simple zero at each of these n points. The case of two or more of these points becoming 
coincident can be treated by taking a limit. The helds cf) and a are smooth functions of 
space and time, so the original expression (4.4) for T has no singularities. However, h 
and X are ill-dehned at the (moving) vortex locations, and V/i, Vy, h and x all diverge 
as a vortex is approached. Let us therefore dehne the region of integration E to be that 
obtained from by removing n discs of small radius e, centred at the vortex locations 
x’’(t). This produces an error of order in T. At the end of the calculation, let e ^ 0. 
During the calculation we can neglect any terms in T (but not in the component helds) 
that are 0(e) or smaller. 

Let C'^{e) denote the circular boundary of the rth disc. The terms involving fi and /2 
reduce to line integrals along C'’'(e). There is no contribution from inhnity as fi and /2 
decay exponentially fast. The integral of ^(Vh ■ Vx) can be expressed in terms of a total 
time derivative. We have 

-I / V/i ■ Vx / |-(Vh ■ Vx) d^x-f2 [ (V/I ■ Vx) i" X dl , (4.16) 

taking into account that the boundary C^{e) is moving with velocity x^. The line element 
dl is that along C^{e), and x^ x dl is a scalar (in two dimensions). If we drop the total 
time derivative from T, which does not affect the dynamics, the integral of ^(V/i ■ Vx) is 
simply replaced by integrals over the circles C‘^{e). 

Combining these observations we hnd 

n Tt n TL n 

T=-^ xd^x+^-T (V/i-Vx)x^xdl -yV / f ■ dl . (4.17) 

^ r=ldC-{e) 

To calculate the integrals along C'^(e) we need the expansions of h and x around the vortex 
centre x^. Since cf) is linear in the neighbourhood of x^, to hrst approximation, h has the 
expansion [11] 

h = log |x - x^(t)|^ + a’' + (dl^xi - xl{t)) + (32{x2 - xl{t)) + ... (4.18) 

where , /?[ and /?2 depend on the locations of the other vortices, and vary smoothly 
with time as the vortices move. The contours of |(/)| are approximately circles near x^. 
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and {(31, measures the extent to which the centres of these circles differ from x’’ as \4>\ 
increases, because of the other vortices. 

From (4.18), we hud that near x’’ 


Vh = 


f2{xi-xl) , 2{x2-xI) , 

V |x-x12 ’ |x-x12 


(4.19) 


and 






X — X 


r 12 


X — X' 


(4.20) 

with higher order corrections that can be neglected. Let us introduce a polar angle 9^ 
relative to the (moving) vortex at x’', with = 0 in the positive xi direction. Then on 
d7'-(e) 


'2cos^^ 2sin^^ 

Vh = ( -+ I3{ , -+ (3. 


(4.21) 


and 


'2 cos 6'"' nr \ -r 
h = -{ - +(3{ ]xl - 


2sin0" 


3-(3 2 ] X 2 3-<3^ + e($[cos9^ + $2^in9^] . (4.22) 


Next, suppose that the gauge has been chosen so that on and inside the circle C^{e), the 
phase X of the held (f) is exactly linearly dependent on the polar angle 9'^, that is 


x = 9^+r 


(4.23) 


where ^|J^, which will be referred to as the orientation of vortex r, depends only on time. 
Then, on C^{e) 


Vx = 


sin 6*^ 
e 



(4.24) 


and 


ic — - {xi sin 9'^ — ±2 cos 9'^) +■ 
e 


(4.25) 


We shall see that it is not possible globally to set '0’' = 0, although it would be on a short 
time interval. 
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We can now evalnate the line integrals along C'^(e) in the kinetic energy expression (4.17). 
Note that dl = e(— sin 9'^, cos 6'^) and x dl = e{x\ cos 9"^ + ±2 d-9'^- Therefore 


r*27r 




{Vh ■ Vx)x’^ X dl = / {-PI sin 9^ + PI cos 9^) {x{ cos 9^ + sin 9^) d9^ 

Jo 

= n{Plxl-Plxl) , 


(4.26) 


and 


r.27r 


f ■dl = 


IC^{e) 


1 + e{P 2 sin 9^ + P{ cos 9^) ) ( - {±1 sin 9^ — ±2 cos 9^) + ip'" ]d9 


1 

2 


r-271- 


{-P{ sin 9^" + PI cos {x\ cos 9^" + sin 9^") d9' 


Stt 


-27,p'"-—{P^2il-Plxl) , 


SO 


T = -7 / xd^x +2 tix^{V + {(Jlxl - PIX 2 )) 


(4.27) 

(4.28) 


r=l 


We still need to consider the integral of y. This is not the time derivative of the integral 
of X, since x is mnltivalned and its integral over the plane ill-dehned. However, a gange 
transformation replaces X by X + X where x is single-valned, and the integral of x changes 
by the time derivative of the integral of x- So, np to a total time derivative, the integral of 
X is gange invariant. It is convenient, in this integral, to extend the region of integration 
back to the whole plane, x is 0{j) near the vortices, so the contribntion to its integral 
from the discs of radins e is 0(e) and can be neglected. 

The integral cannot be evalnated directly. Instead, consider its integral over a hnite time 
interval 

[ [ X dt . (4.29) 

Jto Jk^ 

Snppose that the initial and hnal conhgnrations are the same, that is, the vortex locations 
are the same, and the helds are in the same gange. Consider hrst the simple case where jnst 
one vortex moves anticlockwise aronnd a loop which does not enclose any other vortices. 
At a point ontside the loop, x varies bnt there is no net change in x between the initial 
and hnal time. (At inhnity this is trne, and by continnity this resnlt extends to any point 
ontside the loop.) At a point inside the loop, on the other hand, x increases by 27r. This can 
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be verified by deforming the motion of the vortex around the point into an anticlockwise 
motion of the point around the vortex. The integral (4.29) is therefore 27r times the area 
of the loop. For a more general closed vortex trajectory, possibly with self-crossings, the 
integral is 27r times the signed area enclosed. 

The result generalizes. Even if the loop encloses other non-moving vortices, the integral 
is 2ti times the area of the loop. If all n vortices move, the integral is 27r times the sum 
of the areas enclosed by the n vortex loops. Finally, we must allow for the vortices to 
exchange locations. Suppose two vortices move anticlockwise along half-loops, such that 
they exchange places. Then the integral is 2ti times the area enclosed by the loop. All 
these results can be checked by using a specific model for y, for example the phase of the 
complex polynomial P{z) = ~ where 2 ; and Zrit) are the complex numbers 

representing a general point in the plane and the trajectory of the rth vortex. Then it is 
easy to calculate the change in y at ^ due to the vortex motion. Unfortunately, this model 
doesn’t quite satisfy our requirements on the phase at infinity or on the circles C^{e), but 
this can be dealt with easily. 

Now observe that the correct value for the integral of x over space and time, as discussed 
above, can be obtained from the local expression 

n 

+^ 1 ^ 2 ) , (4-30) 

r—1 

whose integral over time gives again 27r times the sum of the (signed) areas enclosed by the 
vortex trajectories. So the expression (4.30) is equal to the integral of x over the plane, 
up to an ignorable time derivative. Therefore T can be rewritten as 

i'i’’' + (K + ^ 2 ) 4 ; - W'l + . ( 4 . 31 ) 

and this is our final expression for T. The term (4.30) has occurred before, in the context 
of ungauged vortices [20], although the neat calculation in [20] involves the manipulation 
of a divergent integral. For well separated vortices, the terms in (4.31) dominate, since 
I3l and are exponentially small. The coefficients and P 2 have a significant effect 
when the vortex cores are overlapping. Their appearance is not surprising as they also 
appear in the expression for the Riemannian metric on and Samols has determined 
some of their properties which will not all be used here, but which could be useful [11]. 

The terms are topological. Along a closed path in M^, with the initial and final 
configurations identical, the initial and final values of the orientations of the vortices 
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are geometrically the same, but they may have been permuted, and they may also have 
changed by integer multiples of 27r. The motion of the vortices relative to each other 
dehnes a braid, and the integral over time of X]r=i depend only on this braid. In 

fact, each positive generator of the braid group contributes 27r to the integral, as we shall 
see below. 

There is an important consistency check on the expression (4.31) for T. Note that, as 
anticipated, T dehnes a connection on the moduli space The integral of this connection 
around any closed loop on M'^ is called the holonomy around that loop. It can be calculated 
by integrating T over time, for motion (at any speed) around the loop. Since M'^ has no 
singularities, the holonomy should vanish as a loop contracts to a point. In particular, 
consider the closed loop on where two vortices exchange places by moving anticlockwise 
along semi-circular trajectories. This motion is a generator of the braid group. The 
holonomy should become zero as the radius shrinks to zero. Let the vortex trajectories 
be x^(t) = i?(cost, sint) and x^(t) = —i?(cost, sint), with 0 < t < tt. The complex 
polynomial representation P{z) — {z — zi{t)){z — Z 2 {t)) shows that the orientations of the 
two vortices are '^^(t) = t and '^^(t) = 7r-|-t, so the integral of is 27r. Also, by circular 

and reflection symmetry, (/3j, (3]) is of the form /3(i?)(cosf, sint) and (/3f, /3|) = —(/?i, f^D- 
A simple calculation shows that the total holonomy is 

— 2ti + 2'n(3{R)R) . (4.32) 

The coefficient (3{R) is dehned as that occurring in the expansion of h around x^ = (i?, 0) 
when the vortices are at {R,0) and {—R,0). By reflection symmetry, this expansion has 
the form 

h = log((a:i - R)‘^ + a(-R) + P{R){xi - R) + ... . (4.33) 

For large R, the vortex at {—R, 0) has an exponentially small effect, so I3{R) is exponentially 
small and the holonomy is proportional to nR^ — 27r, the 27r being a topological correction 
to the area of the circle enclosed by the vortex trajectories. When R is small, h can be 
estimated simply from the leading logarithmic terms due to both vortices 

h = log((a:i - A)^ ( 0 : 2 )^) + log((a:i -f i?)^ ( 2 : 2 )^) ... . (4.34) 

Expanding about (i?, 0), this becomes 

h = log((2;i - R)‘^ ( 2 : 2 )^) + log4i?^ ^( 2:1 - i?) -^ ... , (4.35) 

H 
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so (3{R) = for small R. This singular behaviour is just what is needed for the holonomy 
to vanish as R 0. 

Presumably, a more sophisticated version of this argument would establish that the holon¬ 
omy vanishes for any loop on as the loop shrinks to a point. 

5. Vortex Motion 

The Lagrangian derived in section 4 for the motion of n vortices is 


L = 27r7 


n 

E 

r=l 


V + + 2 ^ 2 )*! - 








( 6 . 1 ) 


where V is the integral expression (4.3). The term '0^ is locally a time derivative, and 
so has no effect on the vortex motion, although it has a topological signihcance as we 
have seen. A rather similar Lagrangian has been obtained for well separated vortices by 
Dziarmaga [21], using Berry phase methods. 

The general form of the equation of motion for vortex r is 


d f dL 


dt \ dxl 


dL 

dx^ 


= 0 


( 8 . 2 ) 


For L as above, this becomes 


27r7 


27r7 


x'\ 


io + 


dPl dp. 


+ 


dxf dx 2 

dp^2 , dPt 


dx% 


+ 


dx\ 


x^ + 


x% 


dPl dP{ 
dx% dxo 


dp^ dp 


dxl dx\ 


x% 


S 

2 1 ys 

tL/1 


+ 


dV 

dxl 

dV 

dx\ 


= 0 


= 0 , 


(5.3) 


with summation over s implied. It is straightforward to check that V is constant along a 
trajectory. This implies that if n vortices are initially well separated, they remain so, and 
similarly if they are initially all close together. 

Symmetries imply further conservation laws. Consider the Lagrangian for motion on 
in the general form (4.1). If ^“(X) is the vector held generating a symmetry (i.e. the Lie 
derivatives C^A and C^V vanish), then A,q(X)^“(X) is a constant of the motion. The 
symmetries here are translations and rotations in the plane. The associated conserved 
quantities that follow from the specihc Lagrangian (5.1) are, respectively. 


n n 

+ 2 ^ 1 ), + (5-4) 
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and 


E ((<31 + + (/35 + ^ 12 ) 4 ) . (5.5) 

Samols [11] has shown that for any locations of the vortices, 

n n 

E'9r = E'55 = '>' (5-s) 

r=^l r=l 

so the conservation of the qnantities (5.4) implies that the centre of the n-vortex system, 
^Er=lX^ is hxed. 

If there is jnst one vortex, it does not move. This is becanse (3i and /32 vanish, by circnlar 
symmetry or by (5.6), and V is independent of the vortex position. 

Two vortices move aronnd each other in a circnlar motion, rather like two flnid vortices of 
eqnal strength. To see this, change coordinates, writing 

+ Ricos 9, sin 9) 

(5.7) 

— i?(cos 0, sin 9) . 

As was remarked in the previons section, (/di,/?!) = “(/^ij/^i) — /9(i?)(cos6*, sin0). The 
potential V depends only on R. In these coordinates, the Lagrangian simplihes to 

L = 27r7 (29 + - R{R + 2^{R))9^ - V{R) . (5.8) 

The term 29 eqnals 'tp^ +'tp‘^. It has been inclnded, althongh it does not affect the eqnations 
of motion, which are 

XO = ^2° = 0 , R = 0 

2.,[±(R^^2Rm))y=y^. 

The centre X^ does not move, and the relative motion is at constant angnlar velocity on 
a hxed circle. 

At large separation, and (3 are exponentially small, so 9 is exponentially small. The 
maximnm angnlar velocity occnrs at a hnite separation of order the vortex core size. At 
small separation, ^ ~ 0{R^). This is becanse V has a maximnm at i? = 0, and varies 
qnadratically with the “good” radial coordinate on M^, which is R^. The singnlarity 
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in (3{R) does not prodnce a singnlar coefficient of so the angnlar velocity vanishes as 
R tends to zero. In the limit, the two vortices are coincident and at rest. It wonld be a 
nsefnl consistency check to see if there is a corresponding exact, static solntion of the field 
eqnations. 

Symmetry implies that if there are p vortices at the vertices of a regnlar p-gon, and q 
vortices at the centre (possibly with q = 0), then the p-gon will rigidly rotate with its 
centre fixed. 

6 . Conclusions 

In this paper, a time-dependent Ginzbnrg-Landan model for a complex scalar field con- 
pled to electromagnetism in two space dimensions has been considered. The Lagrangian 
incorporates the standard potential energy of Ginzbnrg-Landan theory, and there is a 
Schrodinger kinetic term for the scalar field and a Ghern-Simons term for the electromag¬ 
netic field. The model is exactly Galilean invariant, so fields respond to a transport cnrrent 
simply by a velocity boost parallel to the cnrrent. The model has vortex solntions, and 
at the critical conpling A = 1, it has the 2n-dimensional manifold M'^ of static n-vortex 
solntions satisfying the Bogomolny eqnations. The parameters of these solntions are the 
positions of the n zeros of the scalar field, which can be identified with the vortex positions. 

If A 7 ^ 1, and there are two or more vortices, then they generally move, bnt withont dissi¬ 
pation. For A close to 1, the field dynamics describing vortex motion can be approximated 
by considering solntions of the Bogomolny eqnations with time-varying parameters. This 
leads to a rednced Lagrangian for motion on with a kinetic term linear in the vortex 
velocities, and a potential term. Remarkably, the kinetic term depends only on local data 
associated with each vortex. It defines a gange potential on which depends smoothly 
on the vortex positions, even as the vortices become coincident. 

For two vortices, the gradient of the potential energy is along the line joining them, bnt 
the motion is at right angles to this, and they orbit each other at constant separation. It 
wonld be interesting to stndy the motion of more than two vortices in this model. This 
will reqnire more detailed compntations of the potential energy and the gange potential 
on or a n um erical simnlation of the field dynamics. 

The qnantisation of the vortex motion shonld also be considered. 
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